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Abstract 

Primordial black holes (PBHs) accumulate weakly interacting massive particles 
(WIMPs) around them and form ultracompact minihalos (UCMHs), if the WIMP is 
a dominant component of the dark matter (DM). In this paper, we discuss that the 
UCMHs seeded by the PBHs with sub-earth mass enhance the WIMP annihilation 
in the present Universe and can successfully explain the positron and/or electron 
excess in cosmic ray observed by PAMELA/Fermi experiments. The signal is very 
similar to that from a decaying dark matter, which can explain the PAMELA and/or 
Fermi anomaly without conflict with any constraints as long as the decay mode is 
proper. In this scenario, the boost factor can be as large as 10 5 . In addition, we 
discuss testability of our scenario by gamma-ray point source and gravitational- wave 
experiments. 



1 Introduction 



The origin of the dark matter (DM) is one of the most challenging problems in both parti- 
cle physics and cosmology. In the context of beyond the Standard Model (BSM), a weakly 
interacting massive particle (WIMP) is a good candidate of the DM. In many models, DM 
is considered as thermal relic, whose abundance is determined by its annihilation cross 
section: 



where (av) is thermal average of a product of the DM annihilation cross section and 
velocity. To explain the present DM abundance, Qj^h 2 ~ 0.1 p], the annihilation cross 
section is required to be 3 x 10 -26 cm 3 s _1 , which is a reasonable value for mwiMP = 0(100) 
GeV in many BSMs. 

The annihilation processes also occur in the galactic halo in the present Universe and 
can be a source of the cosmic ray. Recently, the PAMELA collaboration has reported the 
excess of the cosmic ray positron [2] . If the origin of the anomaly is the DM annihilation, 
the required annihilation cross section is very large, typically (av) ~ 10~ 24 cm 3 s _1 , which 
is (9(100 — 1000) times larger than the cross section that accounts for the present DM 
abundance as seen in Eq. ([I]). Therefore, one needs to enhance the annihilation processes 
to explain the PAMELA anomaly. This enhancement is parametrized by so-called "boost 



There are two kinds of the boost factor. One is originated from particle physics. For 
example, the large cross section can be achieved through the Sommerfeld enhancement 
[3] or the Breight-Wigner enhancement [I]. In such cases, the annihilation cross section 
has large velocity dependence and the large cross section is realized at lower velocities. 
Another possibility in particle physics is non-thermal production of the DM. In this case, 
the DM is produced by non-thermal processes such as decay of heavy particles. Therefore, 
the annihilation cross section does not determine the DM abundance as in Eq. ([T]) and 
can have a large value. However, there are severe constraints in both cases since a large 
amount of energy injected by the DM annihilation with such a large cross section can 
change the history of the Universe, such as the Big-Bang Nucleosynthesis (BBN) [5] or 




(1) 



factor. 
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Cosmic Microwave Background (CMB) [6]. In addition, gamma ray observations give 
strong constraints on such a DM [TJ [8] . 

The other source of the boost factor is originated from astrophysics. The annihilation 
rate is proportional to np M , where «dm is number density of the DM. Therefore some 
clumps, in which the DM density is higher than one of the galactic halo, if any, can be 
a source of the boost factor. A candidate of such clumps is DM subhalos, which are 
predicted to be contained in the galactic halo by numerical simulations (9j [TO] • A large 
subhalo, however, has to lie near the Solar System to explain the PAMELA anomaly and 
such cases are improbable [11] . 

In this paper, we propose a new source of the boost factor, which is originated from 
minisize halos seeded by primordial black holes (PBHs) [12j [13], called ultracompact 
minihalos (UCMHs) [HI [151 EE] • The PBHs can be formed in the very early Universe if 
there exist density fluctuations of order of unity. If they make up only a small amount 
of the DM at the formation time, they subsequently accrete the surrounding WIMP DM 
and form minihalos, UCMHs, around them. The WIMP density in the UCMHs is so high 
that a small fraction of them can lead to significant enhancement of the cosmic ray signal 
[TTt [T8] . The boost factor is determined by the UCMH abundance and can be as large as 
O(10 5 ) if there are the UCMHs whose abundance is comparable to that of the DM in the 
present Universe. In this scenario, total flux of the cosmic ray is given by integration of 
each UCMH's flux. Hence the signal is very similar to that from a decaying DM, which 
has no conflict with the BBN, CMB [19] and the gamma ray observations [201 12T] . 

This paper is organized as follows. In Section [2] we describe the annihilation enhance- 
ment in the UCMHs and estimate the annihilation rate per UCMH. In Section [3] we 
estimate the boost factor induced by the UCMHs and give the UCMH abundance nec- 
essary to explain the PAMELA anomaly. Section H] discusses prospects and constraints 
from gamma-ray point source observations. Our summary and conclusion are presented 
in Section [5j 
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2 Annihilation Enhancement in UCMHs 



First, we introduce a WIMP density profile in the UCMH seeded by the PBH and estimate 
the annihilation rate per UCMH, T, based on the discussion in Ref. |17j . 

2.1 WIMP Density Profile in the UCMH 

The PBHs formed in the radiation-dominated era have masses of order of the horizon 
mass at their formation time: 

where T is temperature of the Universal The PBHs lighter than 10 15 g have evaporated 
away by now due to the Hawking radiation, hence, they have to be heavier than 10 15 g 
to remain in the present Universe. The PBH accumulates the WIMP DM on the order 
of its original mass until the matter-radiation equality, z = z eq , and the accreted mass 
increases as oc (1 + z)~ l during the matter-dominated era [HJ EE 122] : 

M h (z) = Mpbh (^^-) • (3) 

After the structure formation, z ~ 30, the growth of the UCMHs proceeds depending on 
the environment around them. The accretion can also stop at l + z$ = /pBH,i(l + £eq)/(l — 
/pbh,i)) when almost all of the WIMPs are accreted. Here, /pbh.i is the initial fraction of 
the PBHs in the DM energy density, which should be much less than unity to accrete the 
WIMPs initially. In this case, the UCMHs constitute the present DM. Note that the DM 
annihilation cross section should be of order of the thermal one even in this case because 
the UCMH mass is dominated by the WIMPs. 

The halo radius at the redshift z, where the WIMP density profile truncates, is given 
by M US 122] 

M h (z)\ l/3 f l + z 



«-w = UAU ^ {—J ■ (4) 

We parametrize the WIMP energy density profile in the UCMH as, 

~ Vi (r< B ) 



1 Here, M ffi is the earth mass, 6 x 10 27 g. We also denote the solar mass, 2 x 10 33 g, as M c 
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where we have denoted i?tr(^eq) as R cq , above which the accreted mass exceeds the PBH 
mass. Here, the indices are assumed to satisfy < Pi < 1.5 and v > 1.5. In the inner part 
of the halo, a large part of the halo is made up of the WIMPs accreted at high redshift 
because the WIMPs accreted later pass through there so quickly. Therefore, the orbits of 
the WIMPs in the inner region, r < R eq , where the PBH dominates the mass, are insen- 
sitive to the density profile there. If the angular momentum of the WIMPs is sufficiently 
small, they free-fall into the PBH and the density profile scales as Pwimp oc r _L5 . On 
the other hand, in the outer part, r> R eq , where the halo mass exceeds the PBH mass, 
the WIMPs feel force from the mass accreted at earlier time. In this region, the profile 
is expected to scale as the self-similar one, Pwimp oc r~ 2 - 25 [22]. Hence, our assumption 
is satisfied in this case. The dark matter density at R eq is determined irrespective of the 
values of the PBH mass as p cq = 1.8 x 10 7 GeV/cm 3 , which is O(10 7 ~ s ) larger than the 
local DM density, p & ~ 0.3 GeV/cm 3 . 

2.2 Annihilation Rate 

The WIMP annihilation rate per UCMH, T, is given by 

I rRtr(z s ) 

T = — / 27rr 2 p 2 VIMP (r)(cn;)dr, ( 6 ) 

m WIMP JR m in 

where z s denotes the redshift when the UCMHs stop growing, which we assume to be 
max(3O,Z0), and -R m i n is a minimum-radius cutoff. For the WIMP density profile ([5]), T 
is determined by the DM annihilation rate in the vicinity of r ~ R eq . Therefore, it is 
convenient to rewrite T as, 

r-ixio s W^^^vJ [-j^j ^ 3xl0 - 2 6 cm3s -iJ' ( 7 ) 

where /profile is defined as 

/profit/ f 2 (^)df (r = r/R eq ). (8) 

The profile-dependent integral, / pro fiie, can be determined irrespective of Rt r ( z s) and R m \ n 
except in the case z/j ~ 1.5 as 



8 x 10 9 GeV/cm 3 f ^^-) ( ^Ur^ , (11) 

1 V100 GeV/ \3 x lO" 26 cm 3 / S y V 1 



In the case i/* ~ 1.5, /profile depends on -R mrn only logarithmically; 

/profile * ^ . (10) 

\ -'l-min / 

Here, the minimum-radius cutoff, R m i n , can be estimated by considering a flattening 
of the density profile due to the WIMP annihilation; the annihilation rate is high enough 
to reduce the density in the inner high-dense region [23] • The shallower profile develops 
from a radius where the density becomes comparable to 

raWIMP _ w 1n9 n ^ r/nm 3 ( m WIMP \ [ ((TV) 

(av)(t -ti 

where to — U is time elapsed from the formation of the UCMH, 4 x 10 17 s. Therefore, the 
minimal-radius cutoff, R m i n , can be estimated to be 

2 2 

^ = (PA~* = o 017 ( mwiMP ^ ^ ( ^ ^1 ~ 3 (12) 

R eq \/W ' UOOGeW V3 x 10- 26 cm 3 / S ; ' 1 ; 

which corresponds to /profile = 0(1). 

We should also consider propagation of the electron/positron in the UCMH; they lose 
their energy by scattering the baryons. For this purpose, we compare the radiation length 
of the high energy electron in hydrogen, Xq ~ 60 g/cm 2 [23], with the corresponding 
length of the baryons within the UCMH, 

fB / PWIMpdr ~ 2f B p e qReq , (13) 

for z/j ~ 1.5. Here, we have assumed that the baryons distribute as ps = /bPwimp in the 
vicinity of r ~ R eq . Substituting Eq. (fl2|) . we obtain 

H min PwrMPdr ^0.13 g/cm (^Gev) ( 3 x 10 -26 cm 3 /s ) ^ 

which is smaller than the radiation length, Xq. Hence, we can neglect the energy loss of 
the electron/positron in the UCMH. 

In conclusion, the order of magnitude of T is, 

r - o(w 29 ) s- 1 ( mwiMP Y 2 ( ( {av) \ (i5) 

1 - U{W ) s ^ ) ^ Mffi x 1Q _ 26 cm3s ^ ) • (15) 
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In the next section, we estimate the abundance of the UCMHs required to explain the 
PAMELA anomaly by using this rate. 

Before ending this section, we have to comment on feasibility of our setup. Though 
we have considered only the UCMHs seeded by the PBHs as sources of the cosmic ray, 
a number of the UCMHs could be formed without producing PBHs [IB] . Density fluc- 
tuations with amplitude > O(10~ 3 ) grow in the matter-dominated era and can collapse 
before z ~ 10 3 . These fluctuations are too small to produce the PBHs, but can seed a 
much larger number of the UCMHs than the PBHs do. Because these UCMHs have no 
core, the self-similar density profile, Pwimp oc r~ 2 25 , is expected to extend to the inner 
region. In this case, the annihilation rate is estimated to be T ~ Pa-^mm- The minimum- 

_ 4 

radius cutoff scales as -R m i n oc p a 9 , so that the annihilation rate strongly varies in time 
as T oc (1 + z). Then, the rate is severely constrained by the recent Fermi observation of 
the extra galactic gamma ray and would moreover violate the constraints from the CMB 
observations [T9] . 

This problem can be avoided if the PBH mass is so small that the accompanying 
small-scale fluctuations which would grow to form the coreless UCMHs eventually are 
wiped out due to collisional damping and free streaming of the WIMPs [2S|. Though the 
collisional damping and the free streaming have comparable length scales, the latter is 
larger by a factor of O(10). The comoving length scale of the free streaming is given by 

/ m WIMP \-| / T kd \ 2 
/fs " 1 PC llOOGeV ) (jOMevJ ' ^ 

where Tkd is the temperature at kinetic decoupling of the WIMPs. The value of Tkd 

depends on the model, but is typically Tkd ~ 10 MeV. On the other hand, a PBH with 

mass Mpbh is formed by the gravitational collapse of the density fluctuations with a 

comoving scale, 

Zpbh * 2 x 10- 3 pc 2 • (17) 

Therefore, if the PBH mass is sufficiently small, 

MmH<2xW ^^y^y\ (18) 
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possible constraints from the CMB observations can be evaded. 

We should also check whether the angular momentum of the WIMPs can be neglected 
because the equations in this section have been derived in the approximation of quasi- 
radial infall. The WIMPs have velocity dispersions induced gravitationally and thermally. 
The former is given by, 

°>^S>,- l mik > (19) 

at comoving scale I where "H is the comoving Hubble parameter and P(k) is the matter 
power spectrum [26]. Hence, it is determined by the matter power spectrum on scales 
smaller than I. The WIMPs turned around at comoving radius, 

^3< 1 + ^0, P c(^)-*(^ (20) 

In the relevant region r ~ R eq , a large part of the halo is made up of the WIMPs accreted 
at z ~ z cq . Therefore, if the PBH mass is taken as, 

(21) 

the gravitational velocity dispersion is exponentially suppressed due to the free streaming 
of the WIMPs and can be safely neglected. 

Next, we consider the thermal velocity dispersion. The WIMPs have the thermal 
velocity dispersion A/3T kd /mwiMP at the kinetic decoupling and the velocity varies as 
oc (1 + z) after that. Hence, the thermal velocity dispersion is given by, 

^ 35 Cm/S ( lOOGev) (lOMevJ (jT^J ' ^ 
The orbit can be considered as quasi-radial at r ~ R eq if R eq is larger than the radius of 
the circular motion (3R tT a t ) 2 /GMpbh, hence the PBH mass should satisfy 

3 

Mp B H > 6 x l(T 4 M e ( ^IMM - 1 ( _^L_] ~~ 2 . (23 ) 
PBH ~ ®Vl00GeV/ ViOMeVy 1 ; 

In conclusion, the angular momentum of the WIMPs can be safely neglected if the PBH 

mass lies in the range 10" 4 M e - M e . Note that the WIMP with the thermal velocity 

cr t has pericenter (3-R tr o"t) 2 /2GMpBH, hence they are not absorbed by the central PBH if 

Mpbh < 10 9 M e (m W iMp/100 GeV)- 3 / 4 (T kd /10 MeV)" 3 / 4 . 



3 Electron/Positron from WIMP Annihilation in UCMH 

3.1 Propagation 

The electron/positron propagates through the galaxy obeying the following diffusion equa- 
tion, 

A ( b (E)f) + V • (D(E)Vf) + Q(S, E) = 0, (24) 

where / is the e 1 * 1 number density per unit energy, b(E) = bo(E/lGeV) 2 is the rate of 
energy loss, D(E) = D (E/lGeV) 5 is the diffusion coefficient, and Q(x,E) is the e ± 
source term. 

In our setup, the source term has a similar form as in the decaying DM case |20j ; 

Q(x, E) = n UCMH r^^, (25) 

where diV e ± / dE is the e ± spectrum emitted by a WIMP annihilation and nucMH is the 
number density of the UCMHs. 

In typical propagation models, propagation length of the electron before losing its 
most of the energy, d, is estimated to be 0(1) kpc(£'/100 GeV) - ^ 1 . The e ± flux can be 
written by the solution of the diffusion equation f )24|) as F = cf /Air. As in the decaying 
DM case, the resultant flux is insensitive to the uncertainty in the UCMH profile in the 
host galactic halo. Only their local abundance is relevant to estimate the flux. 

It can occur in the present circumstances that there are only a few or even no UCMHs 
within the distance d from the Solar system unlike in the decaying DM case. In this case, 
the signals are rather similar to those from nearby clumps [27] • To check whether this is 
the case, we estimate here the number of the UCMHs within the e propagation length, 
N r< d- Assuming that the UCMHs distribute proportional to the DM abundance in host 
the galactic halo as M h nucMH = (^ucmh/^c)p©, where ^ucmh is the cosmic abundance 
of the UCMHs, the result is 



N r<d K 5 x ^ j ^ j ^ ) , (26) 

which is sufficiently large in the relevant mass range 10 _4 M e — M®. Note that the UCMHs 
have an annihilation rate just necessary to explain the PAMELA anomaly in the nearby 
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clump case, T = (9(10 36 ~ 37 ) s~ x [27], when M PB h ~ 1O 3 ~ 4 M and there exists naturally 
one UCMH within 1 kpc if £\jcmh ~ 2 x lO^ 3 " 2 -* (1+Zg)" 1 . However, the abundance at the 
mass Mpbh ~ 1O 3_4 M is tightly constrained by the CMB observations as ^ucmh ^ 10~ 9 
[28] . Moreover, the UCMHs can be formed without producing PBHs in this mass range 
and they would make the situation worse. Therefore, this low-dense case is improbable. 

In the next subsection, we estimate the UCMH abundance necessary to explain the 
PAMELA anomaly assuming that a large number of the UCMHs are contained within 
the distance d from the Solar system. 

3.2 Positron Excesses induced by the UCMHs 

The analysis can be carried out in just the same way as in the decaying DM case [20], so 
that we do not repeat the detailed analysis. Here, we assume that the UCMHs have the 
monochromatic mass function for simplicity. 

Since there are a large number of the UCMHs within the propagation length, we 
can consider a smooth distribution of the UCMHs. Assuming the UCMH abundance pro- 
portional to the DM abundance as before, the source term (1251) within the e ± propagation 
length is estimated to be 

Q 0i p BH ^£>(1(T 25 ) cm-V^ 

o ( 1 + Zs \ ( m wiMP y 2 f (ay) \ dN e ± (27) 

UCMH V 30 yUOOGeW V3 x lO" 26 cm 3 s-V dE ' 

In the case of DM annihilation without enhancement in the UCMHs, on the other hand, 
the source term is given by 

i-3i mv ,-s„-i ( m wiMP \~ 2 ( (crv) \ dN e ± 



Q n PKfi ~ 1.35 x lO" 31 cm-V 1 " "VJ , 1 (28) 

V0 '5^h V100 GeV/ \3 x lO" 26 cm 3 S - 1 J dE K 1 

Therefore, the boost factor is estimated to be 

BE = O(10 6 ) x fi UCMH (^f) ' (29) 

irrespective of the value of the PBH mass. The boost factor is determined by the UCMH 
abundance. In the relevant range of the PBH mass 10 _4 M e — M e , the UCMHs have mass 
in the range 10 _2 M e — 10 2 M e . In this mass range, the UCMH abundance is constrained 
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by microlensing experiments [291 ttH EQ] as 

[0 DM (l(T 2 M e < M h < lO^Me), 
" UCMH< \0.m DM (10-Me < M h < 10 2 M e ). (30) 

Therefore, the boost factor can be as large as 0(10^) in this scenanog 

The signals from the source (j2"7|) is similar to that from the decaying DM with mass 
2mwiMP- To compare with the results obtained for the decaying DM, it would be conve- 
nient to introduce a corresponding decay rate, T c , defined as, 

p _ ^UCMH p 

p /2mwiMP 

which is estimated to be 



or in terms of the DM life-time, r r = F 



-l 

c > 



. . x ^ MH (i±±) ' ( S ) ( 3xl0 tL 5 - ) " ■ (33) 

As seen in Ref. [SHOD], the decaying DM with the life-time, r ~ 10 27 s(m W iMp/100 GeV)- 1 , 
can explain the PAMELA anomaly in the WIMP mass range 10 2-4 GeV without violat- 
ing gamma-ray constraints, though the details depend on the DM decay modes and the 
galactic halo profile. From these, we can conclude that the UCMH abundance 

„ UCMH . W (l±±) - (^) - , ,34) 

is necessary to explain the PAMELA anomaly. Note that the high-dense condition, 
N r< d ^> 1, is satisfied for these values of the UCMH abundance. 

For the WIMP mass mwimp ^ 10 4 GeV, Eq.( l34l) indicates that a small amount of the 



UCMHs is sufficient to explain the PAMELA anomal}^. Hence, Eq. (IM|) can be written 
in a simpler form as 



2 In the case that the accretion stops before z = 30, we can obtain a larger boost factor, O(10 6-7 ), 
as indicated by a factor (1 + z s )/30 in Eq . (f29| . This is because the UCMH mass is smaller by a factor 
30/(1 + z s ) in this case and larger number density of the UCMHs is required to explain their present 
energy density. 

3 Note that (1 + z,;)/30 is always larger than unity. 
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which is much lower than the bounds on the UCMH /PBH abundance [HI [30] for the most 
part of the WIMP mass range. 

On the other hand, the UCMHs have to constitute a large fraction of the DM, ^ucmh — 
Q c , for the WIMP mass mwiMP ~ 10 4 GeV. In this case, Eq. (lM|) reduces to a condition 
on the redshift when the accretion stops, or the initial fraction of the PBHs in the DM 
energy density: 

1-/PBH,* ~ ( HlO*GevJ \3 x 10-26 cmVi ) ' W 
In this case, the UCMH mass should be in the range 10 _8 M e — 10~ 2 M e , where com- 
pact objects like the UCMHs are still allowed to be the DM p3J [30]. For the WIMP 
mass wwimp ~ 10 4 GeV, this UCMH mass range corresponds to the PBH mass range 
10" 10 M ffi - 10~ 4 M ffi . In a lower part of the mass range, we should consider the angular 
momentum of the WIMPs to obtain a more accurate result. 

4 Prospects and Constraints: Gamma-ray point source 

We enclose our paper with a comment on a possibility to observe or have observed the 
UCMHs as gamma ray point sources; the UCMHs are so bright and sparsely distributed 
that they could be observed as gamma ray point sources. 

The flux of the gamma ray from the UCMH is given by <3> 7 = k^T/Airr 2 , where r 
is distance to the point source. We have also denoted the expected number of photons 
per an annihilation process as fc 7 , which has values in the range 0.1 — 10 2 depending 
on annihilation mode and energy range of the photons. The UCMH can be observed as 
a point source if it lies within d po i nt = (/c 7 r/47r$ obs ) 1 / 2 , where <3> b s is the point source 
sensitivity of gamma-ray observations. Then, expected number of the UCMHs observed 
as point sources is estimated as R = (47rdp 0int /3)nucMH- Assuming the annihilation rate, 
T, to be large enough to explain the PAMELA anomaly, we can estimate this number as 

1 J VlO" 8 cm^s- 1 J V100 GeV/ V / \3 x 10" 26 cm^- 1 J ' 1 ' 
Typical values of the sensitivity, $ bs ; is a few x 10~~ 8 cm -2 s -1 for EGRET [32J and 
a few x 10 _9 cm _2 s _1 for Fermi [33], respectively. Therefore, it is concluded that the 
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UCMHs are unlikely to be observed as point sources by the present experiments. In 
FigJH we have displayed the values of R for leptonic annihilation modes obtained by us- 
ing differential sensitivities of present /future gamma-ray experiments [31]. Though the 
UCMHs is unlikely to be observed by the present experiments as noted above, it will be 
possible by experiments in the future depending on the nature of the WIMPs and the 
PBH mass. 




1(T 

"twiMP [GeV] 




10 r 

TOwimf [GeV] 



Fermi fit 




10" 

mwiMP [GeV] 

Figure 1: The expected number of the UCMHs observed as point sources by present / future 
gamma-ray experiments [31] for leptonic annihilation modes. Solid lines indicate the 
expected number observed by HESS (black) and Fermi (red), and dashed lines the number 
expected for ground-based (black) and space (red) gamma-ray experiments in the future. 
Yellow shaded regions are excluded by the Fermi observations of the isotropic diffuse 



gamma-ray [35J. "e excl." shows the region in which e total flux exceeds the Fermi 
observations [36J. The PBH mass is fixed to be the earth mass, M e . The lower value of 
the PBH mass results in the lower values of R as indicated by Eq.([3] 
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5 Summary and Discussion 



In this paper, we have discussed that a small amount of the PBHs with sub-earth mass 
can account for the large boost factor to explain the PAMELA anomaly. The boost 
factor is determined by the present PBH abundance and can be as large as (9(10 5 ) in our 
scenario if the present PBH abundance is comparable to that of the DM. If the large boost 
factor is originated from particle physics such as the Sommerfeld enhancement, there are 
severe constraints from cosmology. In addition, the recent Fermi observation of the extra 
galactic gamma ray gives strong constraints on such a DM [8J. In contrast, the signal is 
very similar to that from the decaying DM in our scenario, which explain the PAMELA 
anomaly without violating such constraints [21]. We have also discussed prospects and 
constraints from the gamma-ray point source experiments and concluded that the future 
gamma-ray experiments can detect the UCMHs as point sources depending on the nature 
of the WIMPs and the PBH mass. 

In our scenario, the cosmic ray signal from the DM is determined by WIMP nature, 
such as the mass, annihilation cross section and mode, and also the PBH abundance. The 
nature of the WIMP can be clarified by the up-coming collider experiment LHC and ILC 
[371 138] . Moreover, the PBH abundance can be also determined by observation of the 
gravitational wave [3H]- If there exist a fraction of the PBHs with mass 10 _4 M e — M e 
as required in our scenario, the gravitational waves associated with the PBHs can be 
detected by the planned space-based interferometers, such as LISA [30], BBO [H] and 
DECIGO [12] ■ Our scenario is testable by the collider and astrophysical experiments. 
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